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Abstract. We study Neumann functions for divergence form, second order elliptic sys- 
tems with bounded measurable coefficients in a bounded Lipschitz domain or a Lipschitz 
graph domain. We establish existence, uniqueness, and various estimates for the Neumann 
functions under the assumption that weak solutions of the system enjoy interior Holder 
continuity. Also, we establish global pointwise bounds for the Neumann functions under 
the assumption that weak solutions of the system satisfy a certain natural local boundedness 
estimate. Moreover, we prove that such a local boundedness estimate for weak solutions 
of the system is in fact equivalent to the global pointwise bound for the Neumann function. 
We present a unified approach valid for both the scalar and the vectorial cases. 



1. Introduction 

In this article, we are concerned with Neumann functions (or sometimes called Neu- 
mann Green's function) for divergence form, second order elliptic systems with bounded 
measurable coefficients in a bounded Lipschitz domains or a Lipschitz graph domain. More 
precisely, we consider Neumann functions for the mxm elliptic systems 

m m d 

(1.1) Yj Li i u ' '=~YjYi D a (A^(x)D / ,u j ), i = \,...,m 

7=1 ;=1 afi=\ 

in f2, where Q is a bounded Lipschitz domain or a Lipschitz graph domain in W 1 with d > 3. 
Here, we assume that the coefficients are measurable functions defined in the whole space 
R d satisfying the strong ellipticity and the uniform boundedness condition; see Section [2] 
for their precise definitions. We do not assume that the coefficients of the system d 1 - 1 b are 
symmetric. We will later impose some further assumptions on the system ( 11. U in the case 
when m > 1 but not explicitly on its coefficients. 

Analogous to the role of Green's functions in the study of Dirichlet boundary value 
problem of elliptic equations, Neumann functions play a significant role in the study of 
Neumann boundary value problem. By this reason, the Neumann functions are discussed 
in many papers, but however, with only a few exceptions, it is assumed that the coeffi- 
cients and the domains are sufficiently regular. In the case when m = 1, Kenig and Pipher 
|[T6l constructed Neumann functions for the divergence form elliptic equations with L°° 
coefficients and derived various estimates for the Neumann functions in the unit ball B. 
Those estimates are the same sorts of estimates known for the Green's functions as appear 
in IfTTl [171 and are nicely summarized in |[T4l Theorem 1.6.3]. Their methods of proof 
are general enough to allow B to be a bounded star-like Lipschitz domain but however, as 
pointed out in |fl9l , it is not immediately clear whether they also work for general bounded 
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Lipschitz domains. Also, their methods do not seem to work for unbounded domains such 
as the half space. On the other hand, Hofmann and Kim 0121 recently proved existence 
and various (interior) estimates for the Green's function of the system ( 11. lb in arbitrary 
domains under the assumption that weak solutions of the system ( 11.1b satisfy an interior 
Holder continuity estimate. Their result has been complemented by a very recent article 
by Kang and Kim ff3l . where global estimates of Green's functions for the system ( II . lb 
are established under some other (but similar) assumptions. In the case when m = 1, 
the De Giorgi-Moser-Nash theory for weak solutions implies such estimates and thus, in 
particular, they were able to reproduce the related classical results of IfTTI IT71 . 

The goal of this article is to present a unified approach for the construction and esti- 
mates of Neumann functions of the elliptic systems ( 11.11 ) in a bounded Lipschitz domain 
as well as in an unbounded domain above a Lipschitz graph. As a matter of fact, it is ex- 
actly where the strength of our paper lies. By using our unified method, we reproduce the 
estimates for Neumann functions of scalar equations with L°° coefficients in the unit ball 
presented in lfl6l as well as those for systems with C a coefficients in C domains appear- 
ing in a recent article lfT31l . Recently, there have been some interest in studying boundary 
value problems for divergence form elliptic equations with complex L°° coefficients above 
a Lipschitz graph; see e.g., fl]|2][3]. In this context, it is natural to consider Green's func- 
tions and Neumann functions for elliptic systems with L°° coefficients in a Lipschitz graph 
domain. In fact, properties of Green's function investigated in 1121 . for elliptic equations 
whose coefficients are complex perturbations of real L°° coefficients, were used in HI. 
However, we are not even able to find a literature dealing with Neumann functions in the 
half space for scalar elliptic equations with U° coefficients. As we have already mentioned, 
our method also goes through in that case, and in particular, we derive the estimates for the 
Neumann function of the scalar elliptic equations with L°° coefficients in a Lipschitz graph 
domain that corresponds to the estimates in |[T4l Theorem 1.6.3]. We hope this article may 
serve as a reference for the Neumann functions and their properties, so that it may become 
a useful tool for other authors. 

We shall now describe our main result briefly. Let Q. c R d with d > 3 be a bounded 
Lipschitz domain or a Lipschitz graph domain. We first construct the Neumann function 
of the system ( 11.11 ) under the assumption that its weak solutions are locally Holder contin- 
uous. In doing so, we also derive various interior estimates for the Neumann function; see 
Theorem 1 3 . 1 1 and |5TI We then show that if the system d 1 - 1 b has such a property that weak 
solutions of Neumann problems with nice data are locally bounded and satisfy a certain 
natural estimate (see the conditions ( ILBb and ( ILB'l i in Section[3]and[5]), then its Neumann 
function N(x,y) has the following global pointwise bound; see Theorem 1 3 . 6 1 and 15721 

(1.2) \N(x,y)\ < C\x-y\ 2 - d , VjcjeQ, x * y. 

Conversely, if the Neumann function has the above pointwise bound, then we prove that 
the system should satisfy the aforementioned local boundedness property; see Theorem l3.9l 
and !5.2l An immediate consequence of our results combined with the celebrated De Giorgi- 
Moser-Nash theory would be that the Neumann function of scalar elliptic equations (i.e., 
m — 1) enjoy the pointwise estimate ( 11.21 ) if Q is a bounded Lipschitz domain or a Lips- 
chitz graph domain. Moreover, if the coefficients of the system ( 11.11 ) belong to the VMO 
class and Q is a bounded C 1 domain, then W ' p estimates imply the aforementioned local 
boundedness property and thus, we would have the estimate ( 11.2b in that case too. As a 
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matter of fact, in those cases, we also have 

\N(x,y)-N(x',y)\<C\x-x'f\x-y\ 2 - d ^ if x + y and 2\x - x'\ < \x - y\, 
\N(x,y) - N(x,y')\ < C\y - y'f\x - yf^ if x + y and 2\y - y'\ < \x - y\ 

for some fi e (0, 1]; see Remark [3~8l 

The organization of the paper is as follows. In Section |2j we introduce some notation 
and definitions including weak formulations of Neumann problems and the precise defini- 
tion of Neumann functions of the system dl.lt . In Section[3j we state our main theorems 
including existence and global pointwise estimates for Neumann functions in bounded 
Lipschitz domains, and their proofs are presented in Section[4] Sectionals devoted to the 
study of Neumann function in a Lipschitz graph domain. In the appendix we provide the 
proofs of some technical lemmas. 

Finally, a few remarks are in order. This article is, in spirit, very similar to IT2l [T3ll . 
where corresponding results for Green's functions have been established. However, the 
technical details are very different since Neumann boundary condition is more difficult to 
handle than the Dirichlet condition. For instance, in Ifl2l . the Green's functions are con- 
structed in arbitrary domains but here Neumann functions are constructed only in domains 
with Lipschitz boundary. We do not treat the case d — 2 in our paper. In dimension two, 
the Neumann functions should have logarithmic growth and requires some other methods. 
As a matter of fact, our method breaks down and is not applicable in two dimensional case. 
One way to overcome this difficulty is to utilize so-called Neumann heat kernel of the el- 
liptic operator defined in a Lipschitz cylinder Q. x (0, oo) c R 3 . However, this approach 
requires first establishing a pointwise bound for Neumann heat kernel that is sharp enough 
to be integrable in f-variable; see J6l|7l for the treatment of Green's functions of elliptic 
systems in two dimensional domains. This topic will be discussed elsewhere because Neu- 
mann heat kernel is an interesting subject in its own right. After submission of the first 
version of this paper, Taylor et al. 11201 constructed the Green's function for the mixed 
problem for elliptic systems in two dimensions. 

2. Preliminaries 

2.1. Basic Notation. We mainly follow the notation used in lfl2l[T3ll . Let d > 3 be an 
integer. We recall that a function <p : R d ~ 1 —> R is Lipschitz if there exists a constant 
K < oo such that 

\<p(x') - <p(y')\ < K\x' -y'\, Vx',y' e R^ 1 . 

A bounded domain Q c R d is called a Lipschitz domain if dQ locally is given by the graph 
of a Lipschitz function. A domain Q c R d is called a Lipschitz graph domain if 

Q = {x = (x\ x d ) e R d : x d > <p(x')}, 

where ip : R d ~ l — > R is a Lipschitz function. Throughout the entire article, we let CI be a 
Lipschitz domain or Lipschitz graph domain in R d . 

For p > 1 and k a nonnegative integer, we denote by W k ' p (Q.) the usual Sobolev space. 
When O is a Lipschitz domain, we define the space W'' 2 (Q) as the family of all functions 
u G W 12 (Q) satisfying J gn u — in the sense of trace. We warn the reader that the space 
W l ' 2 (Q.) is different from the space W 2 (B) used in lfl6l . By using Rellich-Kondrachov 
compactness theorem, one can easily show that there is a constant C = C{d, Q) such that 

(2.1) \\u\\ L 2 (n) < C\\Du\\ LHQ) , Vm e W U2 (£l). 
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The space T 1,2 (£2) is defined as the family of all weakly differentiable functions u e 
L K whose weak derivatives are functions in L 2 (£l). The space T 1,2 (Q) is endowed 

with the norm 

IMIyMp) := \\u\\ L 2dnd-2) (n) + \\Du\\ L 2(ay 

For a Lipschitz graph domain Q with Lipschitz constant K, it is easy to show (see Appen- 
dix) that the following Sobolev inequality holds: 

(2.2) Nlzwcn) <C(d,K)\\Du\y { n), V« e T U (Q). 

We denote Q«(x) = CI n B R {x) and 2" R (x) = <3Q n fi/j(x) for any ^ > 0. We abbreviate 
£2 S = £2 s (x) and 2" s = .T«(x) if the point x is well understood in the context. We define 
d x = dist(x, dQ.) = inf{|x - y| : y e d£2). 

2.2. Elliptic systems. Let L be an elliptic operator acting on column vector valued func- 
tions u = (it 1 , . . . , M m ) r defined on a subset of R'', in the following way: 

Lu = -D a (A a/i D^b), 

where we use the usual summation convention over repeated indices a,/3 = l,...,d, and 
A a P = A a P(x) aremxm matrix valued functions defined on the whole space W 1 with entries 
Ay that satisfy the strong ellipticity condition 

m d 

(2.3) A«f(x)^>A\£\ 2 :=Aj]Yj\e a \ 2 , V£ 6 R"' d , VxeR d , 

i=l a=l 

for some constant A > and also the uniform boundedness condition 

(2.4) |aJ(x)^| < M\€\\v\, V^i7€R" ,J , VxeR d , 

for some constant M > 0. Notice that the z'-th component of the column vector Lu coincides 
with LjjiJ in ( II. U . The adjoint operator 'L is defined by 

'Lu = -DadA^Dpu), 

where A a P = (A^f; i.e., 'A" 13 = A pa . 

2.3. Neumann boundary value problem. We denote by ADu ■ n the conormal derivative 
of u associated with the operator L; i.e., ;-th component of ADu ■ n is defined by 

(ADu ■ tif = A^D fi u j n a , 

where n = (« i , . . . , iid) T is the outward unit normal to dQ.. Let £ be an open subset of dtl 
and / E L] oc (Q) m and g e L) oc (Y.) m . We shall say that u E Wj£(Cl) m is a weak solution of 

Lu — f in Q., ADu ■ n — g on 2 

if the following identity holds: 

(2.5) f Af D^Dj- f gV= f /V, V^eCr(nuS)"'. 

Observe that ( 12.51 makes sense if / is a vector-valued measure in Q; see part ii) in the 
definition of Neumann function below. We are mostly interested in the case when 2 = d£l. 
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2.3.1. Neumann problem in a bounded Lipschitz domain. Let Q C 9. d be a bounded Lip- 
schitz domain. Notice that the inequality (12. U implies that H :- W l ' 2 (Q) m becomes a 
Hilbert space with the inner product 

(2.6) {u,v)b ■= | D a u l D a v [ . 

Jn 

If we define the bilinear form associated to the operator L as 

(2.7) fi(«,v):= \ A^DpuiDj, 

Jn 

then by ( 12.31 ) and (|2.41 l. the bilinear form B becomes coercive and bounded on H. Observe 
that by the inequality ( 12. Il l and the Sobolev imbedding theorem, we have 

(2.8) NlL^> ( n> < C\\Du\\ L 2 (a) = C\\u\\ H , Vh e = W'^Q)'". 
Let / e L 2rf/(rf+2) (Q) m and g e L 2 (3Q) m satisfy the "compatibility" condition 

(2.9) f/+ f g = 0. 

Jn Jan 

Then, by the inequality ( 12.81 ) and the trace theorem combined with d2.1| ), we find that 

F(«):= ff-u+f gu 
Jn Jon 

is a bounded linear functional on H. Therefore, the Lax-Milgram theorem implies that 
there exists a unique u in H such that B(u,v) = F(v) for all v e H = W ,2 (fi) m . Observe 
that any function v 6 W 1,2 (f2) m is represented as a sum of a function in H and a constant 
vector in R m as follows. 



Jan I Jat 



v =: v + c. 

an 



Notice that the condition ( 12.91 ) implies F(c) = 0. Then the identity B(u, v) = F(v) yields 

(2.10) f Af Dpu'Dj = [ fV+ [ gV, VveW^COr. 
Jn Jn Jon 

Therefore, we have a unique solution u in H = W h2 (Q.)"< of the Neumann problem 

Lu — f in Q, 

ADu ■ n — g on 3Q. 

provided / E ^/W+^^yn an( j g e L 2 (3Q) m satisfy the compatibility condition d2!9l ). 

2.3.2. Neumann problem in a Lipschitz graph domain. Let Q c R d be a Lipschitz graph 
domain and recall the inequality i2.2i . Similar to the bounded Lipschitz domain case, 
H :- Y l,2 (Q.) m becomes a Hilbert space with the inner product ( 12.61 ). Also, the bilinear 
form B in 4277} is coercive and bounded in F 1 ' 2 (Q) m . For any / e L 2f/ /< d+2 )(Q) m , the 
inequality (12.21 ) implies that 

F(y):= f /V 
Jn 

is a bounded linear functional on 5 /I ' 2 (Q)'". Therefore, the Lax-Milgram theorem implies 
that there exists a unique u in F 1,2 (Q)'" such that B(u,v) = F(v) for all v e 5 /1 ' 2 (Q)'"; i.e., 
we have 

(2.11) f A'*DpU J D a y= [ fV, Vvey u (0) m . 
Jn Jn 
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Hereafter, we shall say that u is a unique solution in Y 1,2 (£l) m of the Neumann problem 

Lu = f in Q., 
ADu ■ n = on dQ 

if u e Y h2 {Q) m and satisfies the identity d2.111 i. 

2.4. Neumann function. In the definitions below, N = N(x,y) will be an m x m matrix 
valued function with measurable entries Njj : Q. x Q. — » R. 

2.4.1. Neumann function in a bounded Lipschitz domain. We say that TV is a Neumann 
function of L in a bounded Lipschitz domain O if it satisfies the following properties: 

i) N{-,y) e W]£(Cl) and N(-,y) e W U (Q \ B r (y)) for all y e Q and r > 0. Moreover, 
Jan ^(''^ = in the sense of trace. 

ii) LN(-,y) = S y l in Q and ADN(-,y) ■ n = -t^I on dQ. for all y e Q in the sense 

(2.12) f A*D^,y)Drf + -r±- f <p k = <p k (y), e C°°(Q) m . 
Jn Jan 

iii) For any / = (/',... , / ffl ) r e C i 0O (Q) m , the function « given by 

(2.13) «(*):= f N{y,x) T f{y)dy 

Jo. 



is a unique solution in W (Q) m of the problem 

'Lu =f in Q, 



(2.14) 



'ADu ■ n = — f f on dil. 



2.4.2. Neumann function in a Lipschitz graph domain. We say that N is a Neumann func- 
tion of L in a Lipschitz graph domain Q. if it satisfies the following properties: 

i) N(-,y) € Wj£(Cl) and N(-,y) e Y U2 (Q. \ B r (y)) for all y e Q. and r > 0. 

ii) LN(-,y) = 6 y I in Q. and ADN(-,y) ■ n — on 3£2 for all y e Q in the sense 

(2.15) J A* D p N jk (;y)Dj = /(y), = (cf> 1 . . . ,0'") r e C~(0) m . 

iii) For any / e C^(Q)'", the function h given by ( I2.13l l is a unique solution in Y l,2 (Q) m 
of the problem 



(2.16) 



'Lu =f in O, 
'ADu ■ n = on 3D. 



We point out that part iii) in the above definitions give the uniqueness of a Neumann 
function. Indeed, let N(x,y) is another function satisfying the above properties. Then by 
the uniqueness, we have 

f(N- N)(y, xff(y) dy = 0, V/ e («)", 
Jn 

and thus we conclude that N — N a.e. in Q x Q. 
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3. Main results 

The following "interior Holder continuity" condition (IH) means that weak solutions of 
Lu — and 'Lu = enjoy interior Holder continuity. In the case m - 1, it is a consequence 
of the celebrated De Giorgi-Moser-Nash theorem. If m > 1 and d > 2, it is not true in 
general, but however, if the coefficients of the system (II. U belong to the class of VMO and 
if Q is bounded, then it is known that the condition (lIHb holds in that case; see e.g., Ifl2l 
Lemma 5.3]. 

Condition (IH). There exist p e (0, 1] and Co > such that for all x e Q. and < R < d x , 
where d x = dist(x, <9£2), the following holds: If u e W 1 ' 2 (Br(x)) is a weak solution of either 
Lu — or 'Lu — in Br — Br(x), then u is Holder continuous in Br with the following 
estimate: 

(IH) Wm(B m) < C R-«> | \u\ 2 

where [u]cm(B R , 2 ) denotes the usual Holder seminorm. 

Theorem 3.1. Let Q c M. d (d > 3) be a bounded Lipschitz domain. Assume the condition 
(IH). Then there exist Neumann functions N(x, y) of L and N(x, y) of L in Q. We have 
N(-,y), N(-,y) € C^° f (Q \ {y})for ally e Q and the following the identity holds: 

(3.2) N(x,y) := N(y,x) T , Vx,yeCl, x + y. 

Moreover, for any f e L«(Q) m with q > d/2 and g € L 2 (<3Q) m satisfying j Q f + f dn g = 0, 
the function u given by 

(3.3) u(x):= f N(x,y)f(y)dy+ f N(x,y)g(y)do-(y) 

Jn Jan 

is a unique solution in W 1,2 (Q) m of the problem 

{Lu - f in Q, 
ADu ■ n = g on oil. 

Furthermore, the following estimates hold for all y € Q: 
i) \W,y)\\ + \\DN(;y)\\ L i mBilv)) < Cr^ 2 for all r e (0,d y ). 

U) \W,y)\k/M)) < Cr 2 - d+d IP for all r 6 (0,^), where p € [1, ^). 

Hi) \{x e CI : \N(x,y)\ > t}\ < Cr d/( - d - 2) for all t > d 2 r d . 

iv) \\DN{;y)\\ mBr{y)) < Cr l - d+d /" for all r 6 (0, d y ), where pe[l, 

v) \{x € Q : \D x N(x,y)\ > t}\ < Cr^ 1 ' for all t > d\r d . 

vi) \N(x,y)\ < C\x -y\ 2 ~ d whenever < \x - y\ < d y /2. 

vii) \N(x,y) - N(x',y)\ < C\x - x'f°\x - y\ 2 - d -f« if x + y and 2\x - x'\ < \x - y\ < dy/2. 

In the above, C — C(d, m, A, M, Q.,po, Co) > and C depends on p as well in ii) and iv). 
The estimates i) - vii) are also valid for N(x, y). 

Remark 3.5. Observe that if / e L q (Q) m , where q > d/2, satisfies f f — 0, then we may 
take g = in ( 13.31 ) and conclude that 

«(*):= f N(x,y)f(y)dy 
Jn 

is a unique solution in W l ' 2 (Q.) m of the problem 

Lu - f in O, 
ADu ■ n = on dQ. 
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Similarly, if g e L 2 (<9Q)"' satisfies J gn g — 0, then we find that 

b(x):= I N(x,y)g(y)dcr(y) 
Jen 

is a unique solution in W l - 2 {Q) m of the problem 

( Lu -Q in Q, 
1 ADm ■ n — g on dQ. 

Also, the following estimates are easy consequences of the identity (13.2b and the estimates 
i) - vii) for N(-, x): 

i) \\N(x, OIIlW-^x*)) + \\DN(x, Olboiw*)) * C ' -1 ^ /2 /° r fl " r e (0. 
||iV(x, OlltfOWx)) < Cr 2 - d+d/ P for all r 6 (0,d x ), where p 6 [1, jL). 
iiij |{y e Q. : [JV(x,y)| > t}\ < Cr dlid - 2) for all t > d\r d . 

iv) \\DN(x, OWiscbAx)) < Cr l - d+d IP for all r e (0, d x ), where p e [1, 

v) \{ye£l: \D y N(x,y)\ > t}\ < Cr dKd - l) for all t > d\r d . 

vi) \N(x,y)\ < C\x-y\ 2 ~ d whenever < \x - y\ < d x /2. 

vii) \N(x,y) - N(x,y')\ < C\y - yT"\x - yl 2 ^ if x + y and 2\y - y'\ <\x-y\< d x /2. 
In particular, we have \N(x,y)\ < C\x — y\ 2 ~ d whenever < \x - y\ < i ma\(d x , d y ). 

The following "local boundedness" condition (ILBb is used to obtain global pointwise 
bounds for the Neumann function N(x,y) of L in SI. Again, in the case m — 1, it is well 
known that the condition ( ILBt holds in bounded Lipschitz domains; see e.g., [18]. In the 
case when m > 1, this condition does not hold in general and requires certain restrictions on 
the coefficients and domains. It can be shown, for example, that if the coefficients belong 
to the VMO class and the domain is bounded and has C 1 boundary, then the condition (ILBb 
holds via W 1 ' 1 ' estimates; see Appendix. 

Condition (LB). There exists a constant C\ > such that the following holds: For any 
/ € C™(Sl) m and g e C°°(dSl) satisfying j Q f + j gn g = 0, let u e W 12 (Q) m be a unique 
weak solution of the problem 

| Lu = f in Q ( 'Lu - f in Q 

\ ADu ■ n = g on dSl 1 'ADu n = g on dSl. 

Then for all x e Q and < R < diam(Q), we have 

(LB) IMlL-dVfltt) !£ Cj (R- dn \\u\\ L 2 {Clli{x)) + R 2 \\f\\ L 

Theorem 3.6. Let SI c M. d (d > 3) be a bounded Lipschitz domain. Assume the condition 
dim and let N(x,y) be the Neumann function of L in SI as constructed in Theorem \3.1\ If 
we further assume the condition dLBt , then we have the following global pointwise bound 
for the Neumann function: 

(3.7) \N(x,y)\ < C\x-yf~ d for all x,yeSl with x + y, 

where C — C(d, m, A, M, SI, C\). Moreover, for ally & SI and < r < diam(Q), we have 

') \\N(;y)\\ L ^-2 Kci \ Br(y)) + \\DN(;y)\\ L 2 (aXBr(y)) < Cr 1 ^ 2 . 

H) \\N(;y)\\ mBr(y)) < Cr 2 - d+d 'P forpe[l, £). 
Hi) |(iefl: \N(x,y)\ > t}\ < Cr d/(d - 2) for all t > 0. 

iv) \\DN(;y)\\ mBAyy) < Cr'- d+d 'P for p e 

v) |(iefl: \D x N(x,y)\ > t}\ < Cr dKd - 1] for all t > 0. 
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In the above, C — C(d, m, A, M, £2, C\) > and C depends on p as well in ii) and iv). The 
estimates i) — v) are also valid for the Neumann function N of the adjoint 'L. 

Remark 3.8. As we have pointed out, the condition (ILBb is satisfied, for example, in the 
scalar case and also in the case when the system has VMO coefficients and the domain is 
of class C l . In fact, in those cases, we also have the following "local Holder continuity" 
condition: There exist constants //o e (0, 1] and C\ > such that for all x e Cl and 
< R < diam(fi), the following holds: Let u e W 1 ' 2 (Q.r{x))"' be a weak solution of either 

Lu — in Qr(x), ADu ■ n — g on Er(x) 

or 'Lu = in £Ir(x), 'ADu ■ n - g on Er(x), 

where g e C O0 (3Q)" 1 , then we have 

(LH) #n«W s/2W) < Ci (R- d/2 \\u\\ 

By using (ILHb and modifying the proof for the estimate vii) in Theorem l3.1l (c.f. the proof 
for ( 13.7b in Section l4!2b . we have the following global version of the estimate vii): 

\N(x,y) - N(x',y)\ < C\x - x'f°\x - y\ 2 ' d ^° if jc * y and 2\x - x'\ <\x- y\, 

where C = C(d, m, A, M, Q.,fio, C\) > 0. The same estimate is also valid for N. 

Finally, the following theorem says that the converse of Theorem 13. 61 is also true, and 
thus that condition (LB) is equivalent to a global bound (13.7b for the Neumann function. 

Theorem 3.9. Let Q. c E. d (d > 3) be a bounded Lipschitz domain. Assume the condition 
( lIHl l and let N(x,y) be the Neumann function ofL in Q. Suppose there exists a constant C2 
such that we have 

(3.10) \N(x,y)\ < C 2 \x-y\ 2 - d , Vx,yeQ, x + y. 

Then the condition (ILBb is satisfied in £1 with C\ — C\{d, m, A, M, Q, Ci)- 

4. Proofs of main theorems 

4.1. Proof of Theorem fXTJ We closely follow the proof of |[12] Theorem 3.1]. Let us fix 
a function O e C™(R d ) such that <S> is supported in fli(0), < <1) < 2, and L O = 1. Let 
y e Q be fixed but arbitrary. For e > 0, we define 

® s (x) = s- d <l>((x-y)/s). 
Let v = V e ,yji be a unique weak solution in W'' 2 (Q) m of the problem (see Section l2~3l 

I Lv — Q> e ek in Q 
(4 ' !) \ADvn = -(l/\dQ.\)e k on dQ, 

where is the k-th unit vector in W". We define the "mollified Neumann function" 

N s (-,y) = (N E jk (;y))J k=l by 

C 4 - 2 ) ■\; i <--v)-v / -< vr 
Then N e (-,y) satisfies the following identity (see ( 12.101 )): 

(4.3) f A°fD N B jk {-,y)D a <[> i + -±- f f=f 4» £ /, € W u (Q) m . 
By the definition of the space W l,2 {Q.), we have in particular the following identity: 

(4.4) f A^DpNli-^D^ 1 = f <X> e /, E W'-W- 

Jn Jn s (y) 
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By taking = v in (14.4b and then using ( 12.3b . Holder's inequality, and ( 12.8b . we get 



r ®y 

^0.60 



< C e (2 -' /)/2 ||v|| LM / W - 2 , (n) < Cs {2 - d)l2 \\Dv\\ L i (sly 



Therefore, we have (recall v is the fc-th column of N B (-,y)) 

(4.5) l|£W E (-,y)llz2 ( n) < Ce (2 - d)l2 , where C = C(d, m, A, M, Q). 

Let R e (0, dy) be arbitrary, but fixed. Assume that / E C™(£l) m is supported in B R = 
B R (y) c Q. Let u be a unique weak solution in W 1,2 (f2) m of the problem ( 12.14b . We then 
have the following identity (recall 'A?" = A"?): 

(4.6) f A^Dpw'Dj = f /V, Vh> E W 1,2 (Q) m . 
Jn Jn 

Then by setting = it in ( |4.4b and setting w = v £jV;J s in (14.61) . we get 

(4.7) f N? k (;y)f= f <X> E w*. 

Jn Jn s (y) 

Also, by taking w — u in d4.6l ). and using (12.3b , ( 12.8b . and Holder's inequality, we get 

MDu\\l Hn) < ^A^Dpu'Dj = JT/V < C||/|| iM /(* 2)(n) ||D M |b (n) . 
Therefore, we have the estimate 

(4.8) \\Du\\ LHn) < C\\f\\ L ^2 KCl) . 

We remark that the condition (IH) is equivalent to the property (H) in lfl2l Defini- 
tion 2.1]; see lfl2l Lemma 2.3 and 2.4]. By utilizing the condition (IH) and ( 14.8b . and 
following literally the same steps used in deriving Ifl2l Eq. (3.15)], we obtain 

II«IIl~( %2 , < CR 2 \\f\\ L ~ (BR) , 

where C = C{d,m,A,M,Q.,/j.Q,Co). Since R E (Q,d y ) is arbitrary, we get from the above 
estimate and ( 14.7b that 



f W>y)f 



< CR 2 \\f\\ L ~ (BR ), V/eC(B 8 ), Vse(0,/?/2), VRe(0,dy). 

Therefore, by duality, we conclude that 

\\N s (;y)\y (BB (y)) < CR 2 , Vee (0,5/2), VR€(0,d y ). 
Now, for any x E CI such that < \x - y\ < d y /2, let us take R := 2\x - y\/3. Notice 
that if s < R/2, then N E (-,y) E W l ' 2 (B R (x)) m2 and satisfies LN E (-,y) = in B R (x). Then 
by following the same line of argument used in deriving lfl2l Eq. (3.19)], for any x,y £ Q 
satisfying < \x - y\ < d y /2, we have 

(4.9) \N e (x,y)\ < C\x-y\ 2 - d , Ve < \x-y\/3. 

Next, fix any r E (0,d y /2) and let v E be the k-\h column of N E (-,y), where k — 1 m 

and < e < r/6. Let 77 be a smooth function on R d satisfying 

(4.10) < 77 < 1, 77 = 1 on R d \ B r (y), 77 = on B r/2 (y), and \Dt]\ < 4/r. 
We set 4> = ifve in ( 14.4b and then use ( 14.9b to obtain 



(4.11) f 77 2 |Dv e | 2 < C f |£>77| 2 |v £ | 2 < Cr 2 f \x - y\ 2i ^ d) dx < Cr 2 -" 

Jn Jn J B r (y)\B r /2(y) 
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Therefore, by (12.8b and ( 14.111 ). we obtain 

l|v e |lz,MW-2) ( ti\B r cy)) ^ \\W>s\\LWv-V(a) S C||D(77V £ )|| L 2 (n) < Cr (2 ~ d)/2 
provided that < e < r/6. On the other hand, if e > r/6, then ( 14.5b implies 
||v £ || L 2 W - 2)(nWv)) < ||v s || LM / W -2) (n) < C||Dv £ || L 2 (n) < Cr i2 - d)l2 . 
By combining the above two estimates, we obtain 

(4.12) \\N%-,y)\\ L i^-2 K ^B im <Cr i2 - d)l2 , Vr € (0, d y /2), Ve >0. 
Notice from (14.11b and ( 14.10b that for < e < r/6, we have 

\\DN E (;y)\\ L 2 {nXBi<Y)) <Cr (2 -^ 2 . 
In the case when e > r/6, we obtain from ( 14.5b that 

IHW e (-,);)|b ( n\ fl , W ) < \\DN E (;y)\\ LHn) < Cs i2 - d) l 2 < Cr i2 - d)l2 . 
By combining the above two inequalities, we obtain 

(4.13) \\DN%,y)\\ L 2 (nXBriv)) <Cr (2 - d)/2 , Vr e (0, d y /2), Ve >0. 

From the the obvious fact that d y /2 and d y are comparable to each other, we find by (14.12b 
and (14.13b that 

(4.14) \\N s {-,y)\y + \\DN e (;y)\\ L 2 {nWy)) < Cr (2 - d)/2 , Vr € (0, </,,), Ve > 0. 
From d4~T4b it follows that (see £[2] pp. 147-148]) 

(4.15) |{x € Q : \N E (x,y)\ > t}\ < Cr dlid ~ 2 \ Vf > d 2 r a \ Ve > 0, 

(4.16) \{x en-.\D x N s (x,y)\> t}\<Cr d/(d - l \ V/xi^, Vs>0. 

It is routine to derive the following strong type estimates from the above weak type esti- 
mates (flTBT l and (l4~T6b (see e.g. HU p. 148]): 

(4.17) \\N e (;y)\\ mBl . (y)) < Cr 2 - d+d/ P, Vre(0,^.), Ve > 0, for 1 < p < d/(d - 2), 

(4.18) \\DN e (;y)\\ LnB Ay)) ^ Cr l - d+dlp , Vr € (Q,d y ), Ve > 0, for 1 < p < d/(d - 1), 

where C = C(d, m, A, M, Q.,po, Co, p). 

From (14.13b . ( 14.17b . and (14.18b . it follows that that there exists a sequence {e^J^Lj tend- 
ing to zero and a function N(-,y) such that N B,, (-,y) N(-,y) weakly in W l - p (B r (y)) 
for 1 < p < d/(d - 1) and all r e (0,d y ) and also that N B "(-,y) A^(-,y) weakly in 
W h2 (Cl \ B r (y)) for all r e (0, J y ); see El p. 159] for the details. Then it is routine to check 
that N(-,y) satisfies the properties i) and ii) in Section l2.4.1l and also the estimates i) - v) 
in the theorem; see Ifl2l Section 4.1]. 

We now turn to pointwise bound for N{x,y). For any x e Q such that < \x —y\ < d y /2, 
set R :- 2\x - y\/3. Notice that (14.14b implies that N(-,y) e W L2 (B R (x)) and satisfies 
LN{-,y) = weakly in Br(x). Then, by lfl2l Lemma 2.4] and the estimate ii) in the 
theorem, we have 

|iV(*,y)| < CR- d \\N(;y)\\ LHBsix)) < CR- d \\N(;y)\\ LKBiR(y}) < CR 2 d < C\x-yt d . 

We have thus shown that the estimate vi) in the theorem holds. Then, it is routine to see that 
the estimate vii) in the theorem follows from the condition (IH) and the above estimate. 

Next, let x e D. \ {y} be fixed but arbitrary, and let N B '(-, x) e W U (Q)'" 2 be the mollified 
Neumann function of the adjoint operator 'L in Q., where e' > 0. By setting <f> in (14.4b to 
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be the /-th column of N s (•, x) and utilizing an integral identity for N s (•, x) similar to ( 14.4b , 
we obtain the following identity: 



Let N(-, x) be a Neumann function of 'L in Q. obtained by a sequence {sj,}™ j tending to 0. 
Then, by following the same steps as in Ifl2l p. 151], we conclude 

N lk (x,y)=N kl (y,x), V*,J = 1, . . „m, 

which obviously implies the identity ( 13.2b . In fact, by following a similar line of reasoning 
as in [[12, p. 151], we find 

(4.19) N s (x,y) = s- d J ®(iZl ) JN(x,z)dz, 

(4.20) lim N s (x, y) = N(x, y), Vx, y e Q, x * y. 

e—*Q 

Now, let u be a unique solution in W l ' 2 (Q) m of the problem (|2~T4l i with / e C™{Q) m . 
We remark that the condition (IH) implies that u is continuous in Q; see lfT2l Eq. (3.14)]. 
By setting w to be the k-th column of N B (-,y) in ( 14.6b and setting <p = u in (14. 4K we get 



r ^(-,y)/' = r ®y 



We take the limit e — » above to get 

f N ik (x,y)f(x)dx, 



which is equivalent to (12.13b . We have shown that N(x,y) satisfies the property iii) in 
Section l2.4.1l and thus that N(x,y) is a unique Neumann function of the operator L in O. 

Finally, let / e L q (D) m with q > d/2 and g E L 2 (dQ) m satisfy the compatibility con- 
dition (12.9b . and let u be a unique weak solution in W l ' 2 (Q) m of the problem d3.4k see 
Section 12.31 Then u satisfies the identity (12.10b . By setting v to be the fc-th column of 
N e (-, x) in ( 12.10b and utilizing an integral identity for N s (-, x) similar to ( 14.4b , we get 



Jn Jan Jn e (x) 



We remark that the condition ( lIHb together with the assumption that / € L q {Q) m with 
q > d/2 implies that u is Holder continuous in f2; see e.g., Ifl2l Section 3.2]. Then by 
proceeding similarly as above and using ( |3.2b , we obtain 

u\x)= f N ki (x,y)f i (y)dy+ f N^x^gXy) dcr(y), 
Jn Jn 

which is the formula ( 13.3b . The proof is complete. ■ 

4.2. Proof of Theorem 13.61 First, we shall assume that < \x -y| < 2/3 and prove the 
bound ( f3TTb . Let < R < 1 and y e Q be arbitrary, but fixed. Assume that / € C™(Cl) m is 
supported in Qr(j) and let u be a unique weak solution in W'' 2 (Q) m of the problem ( 12.14b . 
Then we have the identities (14. 6b and (14. 7b as in the proof of Theorem 13. II Also, we have 
the estimate ( 14.8b , and thus by ( 12.8b we get 

(4.21) \\u\\ LW -2 Hn) < C\\Du\\i?m < C\\f\\ L ^z m < CR {2+d)/2 \\f\\ L ^n M) , 
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where C = C(d, m, A, M, £2). Observe that g := -(1/|<9Q|) / has the bound 

1*1 1 j^j f Q \f\ ^ CR d \\f\\ L ~ (nMh 
where C = C(d, Q), and thus we have 

(4.22) R\\g\\L^ s( y)) < CR d+1 \\f\\ L ~ (nil(y)) < CR 2 \\f\\ L ~ inM) , 

where we used the assumption that R < 1. Then by (ILBb . ( 14. 21b . ( 14.221 ). and Holder's 
inequality, we obtain 

(4.23) ll"llL»(n s/2W) < CR 2 \\f\\ L - { n M) , 

where C = C(d, m, A, M, Q, Ci). Hence, by d4~7b and ( 14.231 ). we conclude that 



(4.24) 



< C/? 2 ||/|| L » ( n sM) , y/eCM)), Vee(0,*/2). 



Therefore, by duality, we conclude from ( 14.241 ) that 

(4.25) ll# e (-,);)|| L i (£lsW) < CR\ V e e (0,*/2), 

where C = C(d, m, A, M, O, C x ). 

Next, recall that the v = v £ (i.e., £-fh column of N E (-,y)) is a unique weak solution in 
W ia (Q) m of the problem (HIT) . Let x e Q, r > 0, and e > be such that fi £ (y) n B r (*) = 0. 
Then, the condition (1LBI) implies that 

(4.26) II^V £ (-,y)llL»(n,. /2W ) < Ci (r" d/2 ||iV e (-,y)|b ( n rW ) + I^OpV) . 

By a standard iteration argument (see iflOl pp. 80-82]), we then obtain from ( 14.261 ) that 

(4.27) ll# £ (-,>0|| L ~ (a/2W ) < Cr- d W B (;y)\\o ( n r(x)) + Cr, 

where C = C(d,C u \dCl\). 

Now, for any x e Q. satisfying < \x - y\ < 2/3, take R — 3r = 3|x - y|/2. Then by 
( 14.27b and ( 14.25b . we obtain for all e e (0, r) that 

(4.28) |^V £ (x,y)| < Cr- d \\N E (;y)\\ LH[lr(x)) + Cr< Cr- d W t (-,y)\h<p»m + Cr 

< CR 2 ~ d + CR< CR 2 ~ d < C\x - y\ 2 ~ d , 

where C = C(d, m, A, M, Q, C\) and we have again used the assumption that R < 1. There- 
fore, by using ( 14.201 ). we may take the limit s — » in the above inequality and obtain ( 13.7b 
under an extra assumption that \x — y\ < 2/3. In the case when \x — y\ > 2/3, we take 
R = 3r = 1 in d4~28b and get 

\N s (x,y)\ < C < Cdiam(Q) I '- 2 |x-y| 2 ^ < C|x-y| 2 ^, 

where C = C(d, m, A, M, Q, C\). Again, by taking the limit s — > in the above inequality, 
we obtain ( 13.7b even if |x - y\ > 2/3. We have thus shown that (LB) implies ( 13.71 ). 

To derive the estimates i) - v) in the theorem, we need to repeat some steps in the proof 
of Theorem 13. II with a little modification. Let v E be the k-th column of N E (-,y), where 
k = 1, . . . ,m, < s < min(c/ v , r)/6, and < r < diam(£2). Let 77 be a smooth function on 
R d satisfying the conditions ( 14.101 ). We set tp = rfv E in ( 14.3b to get 

(4.29) £ rfA^DpviD^ + jT IrjA^D^D^ + -L jjrj 2 - 1 )v* = 0, 
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-V) 



where we used the fact J gn v e = and // 2 <t> £ = 0. We then use ( 12.3b . d2.4l ), Cauchy's 
inequality to get 

jfV|ZW(-,:y)l 2 < c[fjD V \ 2 \N s (;yf + J-Jjl-^)\N-(.,y)!j, 

where C = C{A, M). By using the conditions in J4. 10b and the pointwise bound for N £ (x,y) 
obtained above, we get 

f \DN B (;y)\ 2 <c(r 2 f \ x - y f-^dx+^- f \x-y\ 2 - d dcr( 

Jci\B r (y) \ JB r (y)\B rl2 (y) J£,-(y) 

< Cr 2 - d , where C = C{d, m, A, M, Q, Ci). 
Therefore, by taking the limit s — > 0, we get 

\\DN(;y)\\ L 2 mBAy)) < Cr (2 - d)/2 , < Vr < diam(Q). 
Observe that the pointwise bound ( 13.71 i together with the above estimate yields 

(4.30) \\N(;y)\\ L ^-2 + \\DN(;y)\\ m a\ Br (y)) < Cr^ d)/2 , < Vr < diam(fi), 

where C = C(d, m, A, M, Q, C\). By following literally the same step used in deriving 
( 14.15b - ( 14.181 ) from ( 14.141 ). and using the fact that \D\ < oo, we obtain the estimates i) - v) 
from ( 14.30b . The proof is complete. ■ 

4.3. Proof of Theorem 13.91 By the symmetry, it is enough to prove (ILBb for weak solu- 
tions of the problem 

( 'Lu=f in Q, 
1 'ADu ■ n — g on dCl, 

where / 6 C?(QJ" and g e C°°(<9Q) m are such that J n f + f gn g = 0. 

Let u be a unique weak solution in W l ' 2 {Q)' n of the above problem. We then have the 
identity (c.f. ( T2~T0l )) 

(4.31) f A"?D/}W j D a u' = f gW + f fw\ Vh> 6 W l ' 2 {Q) m . 
Jn Jan Jn 

Let ( be a smooth function on W 1 satisfying 

(4.32) < £ < 1, supp^ c B R/2 (x), (Elon B 3R/g (x), and \D£\ < 16/R. 
We set w = Cy e in ( 14.311 ). where v e = v e , v ^ is the £-th column of N E (-,y), to get 

(4.33) f {A^DsviDj = - f A^D P UPJ+ f ftV e + f f/V e . 
Jn Jn Jan Jn 

On the other hand, by setting <p = in ( 14.31 ). we get 

(4.34) f *^«*= f ^Af DpviDj ^ f A^DpviD^u' + f 
Jn E oo Jn Jn Jan 

Therefore, by combining ( 14.33b and ( 14.34b , we obtain 

f © s ^= f A^D B viD^u i + -^- f (u k - f A^DbUdJ 
Jnjy) Jn J Jan J n y 



Jan Jn 



NEUMANN FUNCTIONS 



15 



Now, assume that y e Q.r/4(x). Notice from ( 14.32b that dist(y, suppDf) > R/8. Then by 
taking s — > in the above identity, we get 

(4.35) u k (y)= f A^DpNjki-^D^u 1 + — f £u k - f A^D^N jk (;y)D a u i 

Jn Jan Jn 

- f £g i N ik (;y)+ f CfNjki-^-.h+h+h+h + h- 
Jan Jn 

On the other hand, observe that rj 2 v, where v is the k-th column of N(-,y) and r\ satisfies 
the properties in ( 14.101 ). belongs to W 1,2 (£l)'". Then by approximation we may take = rfv 
in ( 12.121 ) to get 

f T?A*D fiV iD a V + f IrjA^D^D^ + f ft 2 - l)v* = 0, 
Jn Jn Jgn 

which corresponds to ( 14.29b in the proof of Theorem l3.6l Following exactly the same steps 
as in the proof of Theorem 13. 61 we then obtain the estimate ( 14.30b . Also, from ( 14.31b and 
the trace theorem, we derive Caccioppoli's inequality 

(4.36) \\Du\\ L 2 ([lm(x)) < C/T 1 HhIIlW)) + CR d,2 (l +R)Mk-<w» + CR d/2+1 \\fh-(n R(x)h 

where C = C(d, m, A, M, Q); see Appendix for the proof. 

Denote A R (y) = Q3 R /4(y) \ B R /g(y). By using Holder's inequality, (14.30b , and (13.101 ). we 
estimate 

\h\ < CR- l \\DN(;y)\\ L 2 (AR(y)) MbdW*)) * CR- d/2 \\u\\ L 2 {nii(x)) , 
\h\ < C\\N(-,y)\\ L 

\I 5 \ < C\\N(-,y)\\ LH a 3m (y))\m\L^n sMx)) < CR 2 \\f\\ L ~ (ClR(JC)) . 
Similarly, by Holder's inequality, the trace theorem, and (14.361 ), we estimate 

\h\ < \dn\- 1/2 Uu\\ L i (m < C(l + /r 1 )||«|b (£M;c)) + C\\Du\\ 

< C(l + R- 1 ) ||K|b (nsW) + C(R d/2 + R d/2+l )\\g\\L^ s (,)) + Cfl' ,/2+1 ||/|| L » ( n fl( *)). 
Also, by Holder's inequality, ( 14.30b . and ( 14.36b . we get 

|/ 3 |< arrive-, y)|| L \\Du\\ L 2 (nm(x)) < CR x - dl2 \\Du\\ L 2 

< CR- d/2 \\u\\ L 2 (nR(x)) +CR(l +R)\\g\h-(w» +^ 2 H/lbxn«(x)). 
Combining together, we get from (14.35b that 

\\u\\ L ~ { n m ) < CR- d/2 (l + R d/2 + R d/2 - l )\\u\\ L 2 (nR) + CR(1 + R dl2 - X + R d/2 + R)\\gh^ R) 
+ CR 2 (l + J?* 2 " 1 )!!/!!^^^)), where C = C(d, m, A, M, Q). 

By a standard covering argument and the fact that R < diam(Q) < oo, we obtain (ILBb from 
the above inequality. The proof is complete. ■ 

5. Neumann functions in Lipschitz graph domain 

This separate section is devoted to the study of Neumann functions in an unbounded 
domain above a Lipschitz graph. 
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5.1. Main results. Since Lipschitz graph domains are necessarily unbounded domains, 
it is more practical to replace the condition (lIHb in Section [3] by the following condition 
dlHl . In the case when the domain is bounded, it is equivalent to the condition dlHt. 
but it is weaker if the domain is unbounded. By the well known De Giorgi-Moser-Nash 
theorem, we have the condition ( lIH't with R c — oo in the scalar case and thus, it reduces to 
the condition i 



Condition (IH'). There exist fi e (0, 1], R c e (0, oo], and C > such that for all x e Q 
and R e (0, d' x ), where d' x : = min^, R c ), the following holds: If u e W 1 - 2 (B R (x)) is a weak 
solution of either Lu — or Lu — in B R — B R (x), then u is Holder continuous in B R with 
the following estimate: 

(IH') Www < C R-x° |h| 2 J . 

Theorem 5.1. Let Q be a Lipschitz graph domain in W 1 (d > 3). Assume the condition 
(lIH'l l. Then there exist Neumann functions N{x,y) of L and N(x,y) of'L in Q satisfying the 
identity (13. 2\ . Furthermore, the estimates i) - vii) in Theorem \3.1\ are valid for N(-,y) and 
N(-,y)for ally e Q provided d x is replaced by d' x — min(c/ c , R c ). 

We also replace the condition ( ILBt in Section |3]by the following condition ( ILB'l i. In 
the scalar case, it is well known that the condition ( ILBl holds in Lipschitz graph domains. 

Condition (LB')- There exists a constant C\ > such that the following holds: For any 
/ e C"(Q) m , let u e F 12 (Q)'" be a unique weak solution of the problem 

Lu — f in £1 ( Lu — f in O 



ADu ■« = onffi 1 <ADu ■ n = on dQ. 

Then for all x e Q and 7? > 0, we have 

(LB') IMIl-w^W) ^ C > (^" d/2 H"llL 2 (n s W) + ^ 2 H/H^(n«(.v))) ■ 

Theorem 5.2. Let Q. be a Lipschitz graph domain in R rf (d > 3) with Lipschitz constant K 
and assume the condition ( lIH'b . If the condition (ILB'l i is also satisfied, then conclusions 
of Theorem 13.61 hold with C — C(d,m, A, M, K,C{). Conversely, suppose there exists a 
constant C2 such that (I3.10l l holds. Then the condition ( ILB'l i is satisfied in Q with C\ — 
Ci{d,m,A,M,K,Ci). 

5.2. Proof of Theorem l5.ll The proof is a slight modification of that of Theorem lXTl Let 
y e Q be fixed but arbitrary. For e > and k = 1, . . . , m, let v = v £ VV t be a unique weak 
solution in Y 1,2 (Q.)'" of the problem 

Lv = E ek in Q, 
ADv ■ n = on dQ, 
and define A^ £ (-,y) by (14.2b . Then A^ E (-,j) satisfies the identity 

(5.4) f A^DpNU-^Dj = f (D £ /, V^F 1 

By the same argument as in the proof of Theorem l3.ll we then obtain ( 14.51 ). Let R e (0, d[) 
be arbitrary, but fixed. Assume that / e C"(Q) m is supported in B R = B R (y) c Q and let u 
be a unique weak solution in T 1,2 (Q)'" of the problem ( 12.161 ). Then, we get the identity d4.7| i 
and also the estimate d4.8l l. By literally the same steps in the proof of Theorem |3.1| we get 
d4.9l l and d4.12| ) - ( 14.201 ) with d y replaced by d' y . Therefore, by the same reasoning as in the 



(5.3) 



' U (Q)" 
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proof of Theorem l3.ll we find that N{-,y) satisfies the properties i) and ii) in Section l2.4.2l 
and also the estimates i) - vii) in Theorem B.ll with d y replaced by d' v ; see lfl2l Section 4. 1]. 
Let u be a unique solution in Y l - 2 (Q) m of the problem (12.16b with / e C~(Q)" ! . Then as in 
the proof of Theorem l3. 1 l again, we get ( 12. 1 3b . Therefore, N(x,y) satisfies the property iii) 
in Section 12.4.21 and thus that N(x,y) is a unique Neumann function of the operator L in 
Q. The proof is complete. ■ 

5.3. Proof of Theorem 15.21 We follow the proofs of Theorem 13.61 and l3~9l with a few 
adjustment. Let y e Q, and R > be arbitrary, but fixed. Assume that / e C™(Q,) m is 
supported in Cluiy) and let u be a unique weak solution in Y 1,2 (Cl) m of the problem 

'Lu = / in Q, 

'ADu ■ n = on <9Q. 

Then we have (14.7b and (14.8b as in the proof of Theorem 13.61 and thus by ( 12.2b we get 
(14.21b . By (ILB'b , (14.21b . and Holder's inequality, we obtain ( 14.23b . Then by following the 
same steps as in the proof of Theorem l3.6l we get ( 14.25b with C = C{d, m, A, M, K, C\). 

Let x e Q, r > 0, and s > be such that B e {y) n B r {x) = 0. Since the A:-th column 
of N s {-,y) is a unique weak solution in F ' 2 (Q) m of the problem (15.3b . the condition (ILB'b 
implies that 

II^V £ (-,y)llL~(n,. /2 (.v,) < C ir - d/2 \\N e (;y)\\ L 2 (nM) . 
Then by following literally the same steps as in the proof of Theorem 13. 61 we obtain the 
desired pointwise bound (13.7b . 

Next, let v £ be the A:-th column of N e (-,y), where k = 1, . . .,m, < s < imn(d', r)/6, 
and r > 0. Let r\ be a smooth function on R d satisfying the conditions in ( 14.10b . We set 
= 77 2 v £ in ( 15.4b to get the Caccioppoli's inequality 

(5.5) f 77 2 |£W e (-,y)| 2 <C f \Dtf\N E (-,y)\ 2 . 
Jn Jn 

By using the conditions in (14.10b for 77 and the pointwise bound (13.7b . we get 

f \DN e (-,y)\ 2 < Cr 1 f |x - y| 4 - 2fl ' c/jc < 2 ~ rf . 

Ja\B,Cv) JB r (y)\B r/2 <y) 

Therefore, by taking the limit e — > 0, we get 

IPtf(-.y)lb<pvM») ^ c/ 2 - d)/2 , Vr > 0. 

The pointwise bound ( 13.7b together with the above estimate yields 

(5.6) \\N(;y)\\ L vm- + \\DN(;y)\\ LHnWy)) < Cr (2 - d)/2 , Vr > 0, 

where C = C(d,m, A,M,K,C\). By following literally the same step used in deriving 
(14.15b - ( 14.18b from ( 14.14b we obtain from ( 15.61 ) the estimates i) - v) in Theorem l3~l6l with 
constants C = C{d, m, A, M, K, C\). 

It remains to show that the pointwise bound (13.7b implies the condition (ILB'b . By the 
symmetry, it is enough to prove ( ILB'b for a weak solution u e Y l - 2 (Q.)'" of the problem 

'Lu =f in 1], 

'ADu ■ n = on dQ.. 

Let u be a unique weak solution in Y x ' 2 (Q) m of the above problem, where / e C™{Q) m , so 
that we have the identity 

(5.7) f AfD^Dj^ f /V, VweF 1 ' 2 (n) m . 
Jn Jn 
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We set w = (v £ in ( 15.7b , where f is as in d4.32l i and v s is the weak solution in Y • 2 (0) m of 
the problem ( 15.31 (i.e., v £ is the A:-th column of N B (-,y)), to get 

(5.8) f tA^Dp^Dj = - f A^Dp^Dj + f f/V,. 
Jn Jn Jn 

On the other hand, by setting = £h in (15.41 ). we get 

(5.9) f <X> e <y= f ^D p v i s D a u i + f A^DpviD^u 1 . 
Jn e (y) Jn Jn 

Therefore, by combining ( 15.81 ) and ( 15.91 ), we obtain 

f O e f h* = f Af DpiD^ii - f A^DpUDaii 1 + f f/V«. 
Jn B (y) Jn Jn Jn 

Assume >' € Q^/4(x) and take s — > in the above identity to get (c.f. ( 14.351 )) 

k*G0 = J A^DpNjki-^D^ij - J" A"jDp^Njk(;y)D 0! u i 

+ f tfNj k (;y)=:h+h + h. 
Jn 

On the other hand, by using the fact that C^(Q) is dense in T 1,2 (£X) (see the proof of 
Lemma 16781 in Appendix), we may set <p - if~v in (12.151 ), where v is the k-th column of 
N(-,y) and rj satisfies the properties in ( 14.10b . to get the following inequality (c.f. (15.51 )): 



f rj 2 \DN(;y)\ 2 <C f |D77| 2 |Af(-,y)| 2 
Jn Jn 



Jn Jn 
Then by proceeding as before, we again obtain the estimate (15.6b . With aid of (12.21 ), we 
also derive the following Caccioppoli's inequality from ( 15.7b : 

(5.10) \\Du\\ L 2 {nm{x)) < CR- l \\u\\ L 2 (ilR(x)) + C\\f\\ L uK^ K n Rm . 

Now, denote A«(j) = Q3«/4(j) \ Br/8(j). Recall that ( satisfies the properties in ( 14.32b . 
Then by Holder's inequality and ( 15.6b . we estimate 

< CR- l \\DN(;y)\\ L 2 (All(y)) ||«lb (ng/jW) < C/T' ,/2 ||h|| lWc)) . 

Similarly, by Holder's inequality, ( 15.6b , and ( 15.101 ). we estimate 

\l2\<CR- l \\N(;y)\\ L 2 \\Du\\ LHnm(x)) < CR- d/2 \\u\\ 

Finally, by Holder's inequality and ( 13.10b . we estimate 

|/ 3 | < C||^V(-,y)|| Ll(n3ff/4tv)) \\fh-(n mix )) < CR 2 \\f\\ L ~ (nit(x)) . 

Combining the above estimates and using a standard covering argument, we obtain ( ILB'I ). 
The proof is complete. ■ 
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6. Appendix 

Lemma 6.1. Let Q. be a bounded C 1 domain in M. d . Suppose that the coefficients A a ? of 
the system ( 11.11 ) belong to the VMO class and satisfy the conditions ( 12.3b and (12.4b . Then 
the condition (ILBb is satisfied. 

Proof Assume that / 6 C~(Q) ffl and g e C°°(dQ)"' satisfy the compatibility condition 
J n f + J an g — and let u e W 1 ' 2 ^)" 1 be a unique weak solution of the problem 

{Lu — f in Q, 
ADh ■ n — g on <9Q. 

Let v = (u, where £ : R d — > R is a smooth function to be chosen later, and observe that v 
is a weak solution in W ,2 (f2) m of the problem 

-D a {A^D p v j ) = £f - ^ - D„n in O, 

A a .^D p vhi a = ft' + F> a on <9Q, 

where we used the notation 

r = A*DrtDpul, F[ y = ,\;fM ; c» ! . 

For i — 1 , . . . , m, let w ! be a solution of the Neumann problem 

( -Aw' = £f - <F in £2, 

{<9w'7<9« = ft'' on dQ. 

Then, by (9j Corollary 9.3] together with the embedding theorems of Sobolev and Besov 
spaces (see e.g., B), we have the following estimate for Dw' provided p > d/(d - 1): 

II-DwIIlccq) < C (||^/Hl^/<p+*(£2) + ll*P|lLP<'/c+<"(n) + llfglb-c-wcan)) • 

Notice that if we set h' a = D a w' + F' a , then v becomes a weak solution of the problem 

f D a (A^D^) = DJi in Q 

| (A^D v j - ti a )n a = on 3Q. 

We then apply (5] Theorem 1] to conclude that v e W 1,p (Q.) m with the estimate 

By choosing £ = 1, we find that « e and 

\\Du\\v(a) < C{\\f\\ D ,dRi»d, {n) + ||^|liPW-w<i(ao)) » 

and thus, via Morrey's imbedding theorem, we find that u e (^(Q) for any fi e (0, 1), 
which particularly implies that u is globally bounded in O. 

To obtain ( ILBb . we employ the standard localization method as follows. Let jefi and 
< R < diam(Q) be arbitrary but fixed. For any yefin Br(x) and < p < r < R, we 
choose the function £ such that 

< ( < 1, supp£ c B r (y), f = 1 on B p (y), and < 2/(r - p). 

Recall that we use the notation 

Q, = Cl r (y) = Q n B r (y), Z r = Z r (y) = dQ. n £ r O0. 
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Then by using the assumptions on £, we estimate terms in (16.21 as follows. 

Kf\\um* m) < Cr 1+dl "\\f\\ L ^n rh 

<C{r-p)- x \\Du\\ 
\\Zg\\w-W(a) < Cr d/ "\\g\\ L ^ r) , 

By using the inequality ( 16.21 ) and the above estimates, we get 

(6.3) \\Du\\ mnp) < 1+I,/ "||/||z.~(n,.) + Cr d/ P\\g\\ L ^ r) 

+ C(r -pr^Klbcn,) + C(r - P r l \\Du\\ 
We fix p > d and let k be the smallest integer such that k > d(l/2 - l/p). We set 

Pi = pd l(d + pi) and r,- = p + (r - p)i/k, i = 0, . . . ,k. 
Then we apply ( 16.3b iteratively to get 

PKlbcn,) < J] C (— j ' ||/|| L „ ( n,) + rf^ \\g\\ L ^ n) ) 

!=I ^ ^' 

+ J] C [— -j ||H|| LP ,- 1( n,) + C* |— -j HDnlIwi^). 
Notice that 1 < p* < 2. By using Holder's inequality we then obtain 

p- d{l/2 - l/p) \\Du\\ LHnp) < Ci -— j (r 1+ ^||/|| L » ( n r) + r d '"\\g\\ L ^ r) ) 

+ c [yr- p ) |, " 1 Nlv W + c(^) / (1/p - 1/2) ||D«ib (£lr) . 

If we take r = /?/4 and p < r/2 = R/4- in the above, then for all y e Qr/^x), we get 

(6.4) p-(Wi-^» |Dh| 2 < CR l+t ""\\f\\ L <°(n K w) + C#<^||g|L~ M 

+ CR-'WuW^n^ + CR d{l/ "- l/2, \\Du\\ L 2 (nmix)) =: A(R). 

Hereafter in the proof, we shall denote Q« = Qr(x). Then by Morrey-Campanato's theo- 
rem (see HU] Section 3.1]), for all z,z' e we have 

\u(z) - u(z')\ < CR [ - d/p A(R), 

where A(R) is as defined in ( 16.4b . Therefore, for any z e Q.r/4 we have 

\u(z)\ < \u{z')\ + \u{z) - u(z')\ < \u(z')\ + CR l - dlp A{R), Vz' e Cl m . 

By taking average over z' £ Qr/4 in the above and using the definition of A(R), we obtain 

sup |«| < -f \u(z')\ dz' + CR 2 \\f\\ L ^ nR) + CR\\g\\ L ~v R) 

+ CR- dl "\\u\\ D , ( n R) + CR x - dl2 \\Du\\ mam) . 
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Then by using Holder's inequality, Caccioppoli's inequality (see Lemma l631 below). and 
the fact that O is bounded, we get 

sup|M| < CR 2 \\f\\ L ~ (nR) + CR\\g\\ L ^ R) + CR- dlp \\u\\ m a R) + CR-' i/2 \\u\\ LHni<) . 

By using a standard argument (see fTOl pp. 80-82]), we derive from the above inequality 
sup|H| < CR 2 \\f\\ L ~ (n , t) + CRMirvu + CR- d/2 \\u\\ L 2 (nii) . 

The proof is complete. ■ 

Lemma 6.5. Let Q c M. d be a bounded Lipschitz domain. Let u e W 1 ' 2 (Q) m be a weak 
solution of the problem 

Lu — f in £1, 
ADu ■ n — g on <3Q, 
where f 6 L°°(Q) m and g 6 L°°(dQ.)" 1 . Then we have 

(6.6) \\Du\\ LHam) < CR- x \\u\\ L 2 (aR) + CR d/2 (l + K)\\gh-(M + CR d/M \\f\k-w 
where C = C(d, m, A, M, Q,). 

Proof. Let rj be a smooth function on R d satisfying 

< 77 < 1, supple fitf, j] = 1 on B R/ 2, and \Drj\ < 4/R. 
By setting v = n 2 u in (12.10b . we obtain 

f rfAfn^Dj^- f InAfDpnu'Dj + f rfgV + f rffu 1 . 
Jn Jn Jen Jn 

Then by ( 12.3b . ( 12.41 ). and Cauchy's inequality, we get 

(6.7) f n 2 \Du\ 2 <C f |£>7 7 | 2 |H| 2 + C||g|| L » (2fi) ||7?H|| Ll( ,n ) + C||/|| L » ( n fl )ll«llL.(n s ). 
Jn Jn 

Observe that the trace theorem (see e.g., (8)) yields 

I \qu\ < C I \D{rju)\ + \tju\ < C \ \Drf\u\ + n\Du\ + \tju\. 
Jan Jn Jn 

Therefore, by using Holder's inequality and Cauchy's inequality, we estimate 
ll*llL-awlMlLi(*i) < Qlgh-to) (^ /2 ^ll«lb(n s ) + CR d,2 \\nDu\y (nR) + C^ /2 ||H|| L2(nR) ) 
< CR d (l + s- 1 + R 2 )\\g\\l^ R) + C7T 2 ||H|| 2 2(n;;) + s f J] 2 \Du\ 2 . 
Similarly, by Holder's inequality and Cauchy's inequality, we obtain 

\\fk-w\M\iHfv < CR d/2 \\fh~(nJ\u\\ LH n R) < CR d+2 \\f\\Un R) + C/T 2 ||H|| 2 2(n;i) . 
By combining ( 16.7b and the above inequality, we get 

||D«|| 2 2(nR/2) < C^- 2 ||«|| 2 2(nfi) + CR d (l + R 2 )\\g\\i^ R) + C^ +2 ||/|||„ (ns) . 
The above inequality obviously yields ( 16.61 ). ■ 

Lemma 6.8. Let fl c R d be a Lipschitz graph domain with Lipschitz constant K. Then, 
for any u e ^'"(Q), we have 

(6.9) IMIzM/Gf-2) (Q) < C(d, K)\\Du\\ L 2 (a) . 
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Proof. We begin with showing that C"(Q) is dense in y 12 (Q). By following the same 
steps as in the proof of approximation theorem for Sobolev functions (see, e.g., (8)), we 
find that C°°(Q) n Y l ' 2 (Q) is dense in Y l ' 2 (£l). On the other hand, one can approximate 
u e C°°(£2) n F 1,2 (0) by a sequence of functions from C£°(Q) in the 7 U (£2) norm as 
follows. For each k = 1, 2, . . ., let e C~(R £/ ) be such that 

< 0* < 1, 0* = 1 on Bfc(O), supp<fo c B 3k (0), and |V0 k | < 1/A:. 

Then, obviously utpt e C^°(Q) and it is easy to check \\wpk - M||yU(n) — > as k — > oo. We 
have thus shown that C^°(0) is dense in F 1,2 (Q). By essentially the same argument, we 
also find that C™(R d ) is dense in F 1,2 (Rr). Therefore, the Sobolev inequality yields that 

(6.10) \\u\\ mM -» m < C(d)\\Du\y m , V M6 y 1 - 2 (R d ). 

Next, we claim that there exists a bounded linear operator E : Y 1,2 (Q) — > 3 /1 ' 2 (R flf ) such 
that Eu — u in Q and 

(6.11) \\D{Eu)\y m < C{d,K)\\Du\\ LHny 

To prove (16.1 11 1, we follow the same steps in the usual proof of extension theorem for 
Sobolev functions in Lipschitz domain (see, e.g., (8)). F° r u e C^°(Q), set 

u + (y) = u(y) if y € Q, 

u~(y) = u(y', 2y(y') - y d ) if y e R d \ Q. 

Note m~ = u + — u on Then, it is routine to check (see (8] Section 4.1]) 

H D "li.W\n) ^ W||D M |b (n) , 

and thus, we have 

ll« _ ll F i.2 (R ^n ) < C(/:)||M|| F i,2 (n) . 

Define 

( u + on Q, 
Eu = u = \ 

{ u on R \ n, 

and note that u is continuous on W 1 . Also, it is easy to see u e Y l,2 (R d ) and 




Du + on Q, 
Dm" on R d \ Q. 



Therefore, we have proved d6.1U in the case when u e C^(£2). Since C"(Q) is dense in 
F 12 (Q), we obtain d6.1U by the standard approximation argument. Finally, we obtain d6.9t 
by combining ( 16.101 ) and ( 16.111 ). ■ 
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